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The measured occupancies of valence orbits in ^^Ge and ^^Se are used as a guideline for modifica- 
tion of the effective mean field energies that results in better description of these quantities. With 
them, in combination with the selfconsitent renormalized quasiparticle random phase approxima- 
tion (SRQRPA) method that ensures conservation of the mean particle number in the correlated 
ground state, we show that the resulting 0u/3f3 nuclear matrix element for the ^^Ge — > ^®Se transition 
is reduced by ~25% compared to the previous QRPA value, and therefore the difference between 
the present approach and the interacting shell model predictions becomes correspondingly smaller. 
Analogous modification of the mean field energies for the A — 82 system also results in a reduction 
of OuPP matrix element for the ®^Se *^Kr transition, making it also closer to the shell model 
prediction. 
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I. INTRODUCTION 

The fundamental importance of the search for Oi^/3/3 
decay is widely accepted (see, e.g. the APS Study of 
Physics of Neutrinos [T]). Observing the decay would 
tell us that the total lepton number is not a conserved 
quantity, and that, consequently, neutrinos are massive 
Majorana fermions. Experimental searches for the Qu(3j3 
decay, of ever increasing sensitivity, are being pursued 
worldwide (for a recent review of the field, see How- 
ever, interpreting existing results as a measurement of the 
neutrino effective mass, and planning new experiments, 
depends crucially on the knowledge of the correspond- 
ing nuclear matrix elements that govern the decay rate. 
Accurate determination of the nuclear matrix elements, 
and a realistic estimate of their uncertainty, is therefore 
an integral part of the study. 

The nuclear matrix elements for Qvfip decay must be 
evaluated using tools of nuclear structure theory. Unfor- 
tunately, there are no observables that could be directly 
linked to the magnitude of Qvj3[3 nuclear matrix elements 
and that could be used to determine them in an essen- 
tially model independent way. In the past, knowledge of 
the 2v (3 P-decay rate and of the ordinary (3 decay ft val- 
ues were used to constrain the nuclear model parameters, 
in particular when the Quasiparticle Random Phase Ap- 
proximation (QRPA) was employed |2I|1]. Clearly, when 
other relevant data become available, and the nuclear 
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model is constrained to reproduce them, confidence in the 
deduced Qv(3j3 nuclear matrix elements would increase. 
Recently a set of such data, the occupation numbers of 
neutron valence orbits in the initial ^^Ge and final ''^Se 
nuclei, were determined in a series of measurements of 
cross sections for neutron adding and removing transfer 
reactions [5] . A similar series of measurements involving 
proton transfer reactions also became recently available 

Here we examine in detail how sensitive the matrix 
elements are to these quantities, and how much the pre- 
viously determined nuclear matrix elements change when 
the input of the nuclear model is modified so that occu- 
pancies of individual orbits are correctly reproduced. As 
in the previous calculations [3] , we use the QRPA method 
and its generalizations. 

The occupation numbers for orbits with angular mo- 
mentum j (and any other quantum numbers) in the ini- 
tial nucleus, measured experimentally, are simply 

and the same quantity is determined for the ground state 
\^%n) °^ ^^"^ ?n\&\ nucleus. Here c~^^ is the creation oper- 
ator for a proton in the orbit jp or a neutron in the orbit 
jn and Cj^m is the corresponding annihilation operator. 
The states \^fnit) ^"^^ lO/m) are the true ground states 
with all correlations in them. 

Theoretically, for pure pairing BCS wave functions, the 
occupation numbers 

- (0^cs|S™4mC,-,„|0+^s) - v] X (2j-+ 1) , (2) 

depend only on the amplitudes Vj or Vj^ that are ob- 
tained by solving the gap equations. The amplitudes v 
are constrained by the requirement that the expectation 
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value of the total neutron and proton numbers are con- 
served, i.e., 



iV(orZ) = E„(p)' 



RCS — 



Mp)'"n(p) 



X (2j„(p) + 1) . (3) 



Here, and in the following, we use n and p to label all 
quantum numbers of the corresponding neutron or pro- 
ton orbits. 

However, in the correlated QRPA ground state the oc- 
cupation numbers are no longer the pure BCS quantities. 
Instead, they depend, in addition, on the solutions of the 
QRPA equations of motion for all multipoles J, and can 
be evaluated using 
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6i(p) = (2in(p) + 1) {OqrpaI 

(5) 

is the expectation value of the number of quasiparticles 
in the orbit n{p). (Here ^ ^ mi'^j„(p)-m £^re the creation 
and annihilitation operators for the quasiparticle with 
quantum numbers n(p),m.) 

The quasiparticle occupation numbers Cn(p) can 
be obtained iteratively using the equations of mo- 
tion of the Renormalized Quasiparticle Random Phase 
Approximation (RQRPA) and of the Selfconsistent 
RQRPA(SRQRPA) through the renormalization factors 

T^pn = 1 - - 6i (6) 

= 1 - 2^^"'^f"' (Sj,fc(2J+ l)|F;„'f ' 



QRPA 



are no 



where = 2?pn' ^^-il"- 

Note that the occupation numbers '^■^(p) 
longer constrained by the same requirement, Eq.(|3|, that 
the particle number is conserved on average, as the BCS 
occupation numbers are. 

The past applications of QRPA to the evaluation 
of the Qvj3f3 nuclear matrix elements used standard 
parametrizations of the nuclear mean field, usually in the 
form of the Coulomb corrected Woods-Saxon potential 
fitted globally to a variety of nuclear properties (typi- 
cally used parameters are those quoted in [T' or in [8 ). 
In few papers [S] attempts were made to modify the sin- 
gle particle energy input in order to better describe the 
energy levels of the /3/3-decay candidate nuclei. It turns 
out that, at least in the A = 76 case, the quantities n^(^J^ 
based on the Woods-Saxon potential single particle en- 
ergies, do not agree well with the experimental results 
of [S] . In particular, the occupancy of the neutron gig /2 
orbit appears to be underestimated. 



Shortly after the data of Ref. [5] became available, 
a new publication |10j appeared, where the neu- 
tron single particle energies (only the valence orbits 
.99/2: /5/21P1/2 E^nd P3/2) were modified so that the neu- 
tron number occupancies of ^^Ge and ''^Se, for which the 
quantities n^^^ were used, were better reproduced. The 
proton mean field energies were also modified so that the 
one-quasiparticle energies in the odd-Z nuclei ''^As and 
^^Br were also better described. The authors of Ref. [TU] 
conclude that this modifications result in sizable reduc- 
tion of the 0z//3/9 nuclear matrix element for the ''^Ge 
^^Se transition. 

In the present work we carefully analyze the role of the 
constraints represented by the knowledge of the orbit oc- 
cupancies n^f^^ . Clearly, these quantities refiect presence 
of correlations beyond pairing correlations described by 
^n(ji) ■ ^ closest substitute for these correlations we 
use here rfif^f"^. Moreover, to describe such correlations, 
we use the SRQRPA that, as the simple BCS, and unlike 
the QRPA or RQRPA, conserves the neutron and proton 
particle numbers on average. (We describe the method 
in more detail in the next section). 

Let us stress that the comparison between the mea- 
sured and calculated occupancies in the experimental 
Refs. [5, .6^, as well as in the theoretical paper [1^, was 
based on equating the experimental values "-^^p) with the 

BCS values "•^(pf • As pointed out above this is not re- 
ally a justified comparison as far as the QRPA and its 
generalizations are concerned. 

The paper is organized as follows. In the next section 
we briefiy describe the SRQRPA method, and illustrate 
the effect of QRPA correlations on the occupation num- 
bers. In Section III we discuss the modifications of the 
mean field energies for the A ~ 76 system that results 
in better description of the occupation numbers of va- 
lence orbits. We also show that using the modified mean 
field energies improves the description of the contribu- 
tion of low-lying states to the 2vf3f3 matrix element. In 
Section IV we present our result for the "^^Ge "^^Se 
Qvf3j3 nuclear matrix element and show that the mean 
field adjustment needed in order to better describe the 
orbit occupancies leads to reduction of the difference be- 
tween the QRPA and nuclear shell model results 
Analogous modifications are applied in Section V to the 
yl = 82 system with two more protons and four more 
neutrons then in A = 76. It also results in a noticeable 
reduction of M^'^ for the ^^Se Qvfip decay. We describe 
the contribution of individual orbits to the Qv(3j3 matrix 
element in Section VI and conclude in Section VII. 



II. SELFCONSISTENT QUASIPARTICLE 
RANDOM PHASE METHOD 



The standard QRPA method consists of two steps. 
First, the like-particle pairing interaction is taken into 
account by employing the quasiparticle representation. 
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In the second step the hnearized equations of motion are 
solved in order to describe small amplitude vibrational- 
like modes around that minimum. In the renormalized 
version of QRPA the violation of the Pauli exclusion prin- 
ciple is partially corrected. 

The drawback of QRPA and RQRPA is the fact that, 
unlike in BCS, and as mentioned already earlier, the par- 
ticle number is not conserved automatically, even on av- 
erage, i.e., 

E^.^n^Q«^-^^iV (7) 

and the same is true for the proton states. Naturally, 
in the limit of negligibly small amplitudes, when the 
quasiparticle occupation numbers £,p and ^„ in Eq. ^ are 
small, the inequality in ([t]) is correspondingly small as 
well. However, for realistic hamiltonians the differences 
between the left-hand and right-hand sides of Eq.(|7| in 
QRPA is of the order of unity (an extra or missing neu- 
tron or proton). 

The selfconsitent renormalized QRPA method 
(SRQRPA) removes this drawback by treating the BCS 
and QRPA vacua simultaneously. For the neutron- 
proton systems, of interest in the present context, the 
method was proposed and tested on the exactly solvable 
simplified models in Refs [T^]. It is a generalization of 
the procedure proposed earlier in [TS]. 

Here we briefly describe the basic features of SRQRPA. 
In QRPA, RQRPA and SRQRPA the phonon operators 
are defined as 

where X^^^^^j and ^(p^y are the usual variational ampli- 
tudes, and ^(p„)jjv/ angular momentum coupled 
two-quasiparticle creation operator. The X and Y am- 
plitudes, as well as the corresponding energy eigenvalues 
Uk are determined by solving the QRPA eigenvalue equa- 
tions for each 

The matrices A and B above are determined by the 
hamiltonian rewritten in terms of the coupled quasipar- 
ticle operators: 

, ,= 

- —1/2 

where ^(p„)j,Af = Dpn A(p„)j,M. The resulting matri- 
ces are independent of the angular momentum projection 
M. 

In RQRPA and SRQRPA the nonvanishing values of 
T^pn — 1 is taken into account by using the amplitudes 

^(pn.J") — -^pn ^(pn,J")' ^ {pn,J^) ~ '^pn ^(pn,J") ' 

(11) 



instead of the standard X and Y, everywhere and also in 
the QRPA equations of motion. 

By doing all of this an inconsistency appears between 
the BCS, with the ground state lO^,^^), and QRPA (as 
well as RQRPA) with the ground state \Oqjipa)- I^i 
SRQRPA this inconsistency is overcome by reformulating 
the BCS equations. This is achieved by recalculating the 
u and V amplitudes. In SRQRPA the state around which 
the vibrational modes occur is no longer the quasiparti- 
cle vacuum, but instead the Bogoliubov transformation 
is chosen is such a way that provides the optimal and 
consistent basis while preserving the form of the phonon 
operator, Eq. ([8|. The modified coefficients of the Bo- 
goliubov transformation still fulfill the basic requirement 
that the so-called dangerous graphs, terms in the Hamil- 
tonian with only two quasiparticle creation or annihila- 
tion operators, vanish. 

In practice, the SRQRPA equations are solved itera- 
tively. One begins with the standard BCS m, v ampli- 
tudes, solves the RQRPA equations of motion and calcu- 
lates the factors Dpn- The u,v amplitudes are recalcu- 
lated and the procedure is repeated until the selfconsis- 
tency is achieved. 

The SRQRPA was applied initially to the evaluation of 
2i//3/3 matrix elements in Ref.fT? and to the evaluation of 
OvfiP matrix elements in Ref.[15j. Numerically, the dou- 
ble iteration procedure represents a challenging problem. 
To simplify it, in Refs. [Ml [15] the bare interaction was 
used, and no attempt was made to fit the odd-even mass 
differences. In addition, no adjustment of the particle- 
particle coupling constant gpp was made, and gpp — 1 
was used. Consequently, the numerical values disagreed 
noticeably with the experiment in the 2vj3[3 case, and 
with calculations by other authors in the Qvf3(3 case. 

The numerical problems were resolved in Ref. |16j 
where instead of the G-matrix based interaction the pair- 
ing part (and only that part) of the problem was replaced 
by a pairing interaction that uses a constant matrix ele- 
ment whose value was adjusted to reproduce the experi- 
mental odd-even mass differences. This is the procedure 
that we adopt also here, after showing that within the 
QRPA the replacement of the G-matrix by a constant 
pairing matrix element makes little difference (see be- 
low). Thus, in the iterative procedure only the chemical 
potentials A„ and \p are changed. 

Adopting this simplification, and using the usual re- 
quirement, as in fj|, namely that the 2i> f3P decay rate is 
correctly reproduced by renormalizing the coupling con- 
stant gpp correspondingly, the authors of Ref. [H] have 
shown that the 0i//?/3 matrix elements evaluated with 
SRQRPA agree quite well with the matrix elements of 
Refs. [3 . The requirement of conserving the particle num- 
ber have not caused substantial changes in the value of 
the Qvj3P matrix elements in that case. 

In Table [l] we illustrate the problem of the particle 
number nonconservation within QRPA and to some ex- 
tent also in RQRPA, and its restoration in SRQRPA. The 
case of ^^Ge ^^Se is chosen, with ''''Ca as a core, and 
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with the p,f and s,d,g shells (9 single particle orbits) 
for both neutrons and protons included. Note that in 
Ref.finj, as noted above, only the BCS occupancies were 
considered. We believe, contrary to the arguments there, 
that the effect of average particle number nonconserva- 
tion in the QRPA vacuum need to be considered. Us- 
ing SRQRPA, which does conserve the average particle 
number, is certainly more consistent when the problem 
of individual orbit occupancies is addressed. 

Even though the conservation of the average particle 
number is almost restored in RQRPA and, as we will 
show below, the numerical values of the M°'' are very 
similar in RQRPA and SRQRPA, we still prefer to use 
the selfconsistent method. Among other things, the vio- 
lation of the Ikeda sum rule, which in the RQRPA is as 
large as 20%, is reduced substantially (but not eliminated 
completely) when SRQRPA is employed. 

TABLE I: The expectation values of the particle number op- 
erator. The mean field energies as in Ref. [3]. For QRPA, 
RQRPA and SRQRPA the particle-particle interaction renor- 
malization constant gpp is chosen from the requirement that 
the 2i^/3/3-decay rate is correctly reproduced. 



System: 


BCS 


QRPA RQRPA SRQRPA 


initial protons 


12.00 


13.05 


12.05 


12.00 


final protons 


14.00 


14.61 


14.01 


14.00 


initial neutrons 


24.00 


23.86 


23.98 


24.00 


final neutrons 


22.00 


22.16 


21.95 


22.00 



Finally, to see the difference in treating the pairing part 
of the problem using the realistic G-matrix based interac- 
tion (but adjusting its strength so that the experimental 
pairing gaps are correctly reproduced, as was done in 
Refs. [3]) and calculation performed with the schematic 
pairing force (a constant matrix element adjusted simi- 
larly) we quote the QRPA and RQRPA values of the M^" 
matrix element, again for the ^^Ge — *■ ''^Se case. With re- 
alistic pairing interaction we obtain M'^'^ = 4.3(3.8) with 
QRPA(RQRPA) while with the schematic pairing inter- 
action the result is 4.4(3.9). Using the schematic pairing 
makes little difference in this case. 

III. THE ROLE OF ORBIT OCCUPANCIES 

The occupancies of the neutron and proton valence or- 
bits in ^^Ge and ^^Se were determined experimentally in 
Refs. [316]. As shown in Section I, within QRPA and its 
generalizations, the occupancies of individual orbits, cor- 
responding to the summed spectroscopic strength mea- 
sured in Refs. [5j [6], are determined not only by the 
BCS amplitudes u and v but also by the quasiparticle 
content of the correlated ground state jOgj^p^). In turn, 
the BCS amplitudes u and v depend on the nuclear mean 
field energies, on pairing gaps A that are fitted to agree 
with the known odd-even mass differences, on the chem- 
ical potentials A that are determined by the requirement 



that the particle number is conserved on average and, 
within SRQRPA that we adopt, indirectly on the solu- 
tions of the phonon equations of motion. 

In Refs. |S] the mean field was based on Coulomb cor- 
rected Woods-Saxon potential using the globally fitted 
parameters quoted in [7]. The resulting valence orbits 
occupancies do not agree very well with experiment in 
that case. In order to describe the experimental occu- 
pancies better, we modify the input mean field energies 
to some extent, mainly for the valence orbits. Since our 
primary goal is to evaluate the nuclear matrix elements 
for the 0i>f3P decay which depend on the quasiparticle 
energies only weakly, modifying the mean field energies 
essentially means that, through the u and v, the occu- 
pancies are adjusted. 

On the other hand, we use the known rate of the 
2i^/3/3-decay to fix the renormalization constant gpp, the 
strength of the neutron-proton particle-particle force. 
The matrix element for the 2v(3P decay, in turn, de- 
pends on the energies of the 1+ states significantly. In 
the intermediate nucleus ^^As the energies of a few low- 
lying 1+ states are known. In calculation that used the 
global Woods-Saxon potential [3^ these energies were not 
described well either, indicating again that the global 
single-particle potential is not optimal for the A = 76 
system. 

Guided by such considerations we modified the mean 
field energy input, determining a set of effective single- 
particle energies for ^^Ge and ''^Se that gives, essentially 
within errors, the measured valence orbit occupancies 
calculated using the SRQRPA. This effective set, at the 
same time, improves the description of the energies of 
low-lying 1+ states in ''^As. In constructing the effec- 
tive single-particle energies, we kept, unlike in Ref. |10j . 
the globally fitted spin-orbit splittings of all orbits in- 
tact. The neutron and proton mean field energies used 
previously in Refs. [3] are compared with the adjusted 
set used further here in Figs, [l] and |2] 

As one could see, the main difference is the overall 
shift of the neutron gds shell closer to the fp shell. The 
proton levels are shifted less, with the gds shell lifted fur- 
ther away from the chemical potential A and from the fp 
shell, i.e., an opposite tendency than in the case of neu- 
trons. The resulting occupancies of the valence subshells 
are shown in Table [IT] and compared with the measured 
values [5] [B] (in columns 5 and 9). The occupancies in 
columns 2 and 6 are those of Refs. ^ evaluated from 
the BCS expression, Eq. with the standard Woods- 
Saxon potential as input. Those are the values quoted in 
the experimental papers [ni[S|- The values in columns 3 
and 7 are the rfi^^^ evaluated with the standard QRPA 
and using the standard Woods-Saxon potential. We dis- 
play them in order to stress the importance of proper 
theoretical treatment, and to show how much difference 
the ground state correlations make. Finally, in columns 
4 and 8 are the occupancies that we use further here and 
which are evaluated in the correlated ground state us- 
ing the SRQRPA and the adjusted set described above. 
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FIG. 1: Comparison of the neutron levels in the Woods-Saxon 
potential used in Refs. [3] (WS) and the adjusted mean field 
energies used here and described in the text. Symbols A in- 
dicate the chemical potential and the crosses indicate the oc- 
cupancy of the individual orbits. 
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FIG. 2: Comparison of the proton levels of the Woods-Saxon 
potential used in [3] and the adjusted mean field energies used 
here and described in the text. Notation as in Fig. [l] 



Note that the entries in columns 3,4,7 and 8 were evalu- 
ated with the correlated QRPA state obtained with the 
coupling constant gpp determined in the usual way, i.e. so 
that the 2i//?/3 rate is correctly reproduced. The sum of 
the corresponding entries, i.e. the total calculated num- 
ber of neutrons and protons above the ""^Ca core were 
already given in Table |T] 

The overall improvement in describing the occupan- 
cies is clearly visible. For neutrons, the total calculated 
occupancy of the valence shells is calculated to be 15.3 
and 13.4 for ^^Ge and ™Se, respectively, with 0.4(0.6) 
neutron vacancies in the fjf2 and 1.0(1.2) neutrons oc- 
cupying the rest of the gds shell. In particular, in the 1^5/2 
there should be 0.35(0.39) neutrons according to our cal- 
culation, while Rcf. 5 tentatively assigns 0.2 neutrons to 



^5/2 in ^^Se. For protons the total occupancy of valence 
shell is calculated to be 4.4 and 6.3, respectively, with 
0.8(0.9) proton vacancies in fj/2 and 0.4(0.6) protons in 
the rest of the gds shell. 



TABLE II: The calculated occupancies of individual neutron 
and proton orbits for the two considered nuclei and using 
BCS only in columns 2 and 6, standard QRPA in columns 
3 and 7 (label Q, these entries were obtained with the stan- 
dard Woods-Saxon potential) and the average nucleon num- 
ber conserving SRQRPA in columns 4 and 8 (label S, these 
entries were obtained with adjusted single particle energies). 
In columns 5 and 9 are the experimental occupancies of va- 
lence orbits [3 [6]. The first entry, for the p orbit, is the sum 
of occupancies of pi/2 and P3/2 orbits. 



neut. 

P 

ihji 

f7/2 
Sl/2 
^3/2 
^5/2 
37/2 
99/2 



prot. 
P 

/5/2 
/7/2 
Sl/2 
^3/2 
^5/2 
37/2 
99/2 



Ge 



BCS 
5.65 
5.54 
7.91 
0.01 
0.03 
0.09 
0.14 
4.63 



Q 
5.27 
5.12 
7.67 
0.05 
0.14 
0.30 
0.53 
4.78 



S 
4.64 
4.34 
7.62 
0.07 
0.15 
0.36 
0.48 
6.35 



exp 
4.9±0.2 
4.6±0.4 



.5±0.3 



2.23 
1.61 
7.83 
0.00 
0.01 
0.01 
0.02 



2.34 
2.27 
7.19 
0.02 
0.07 
0.12 
0.19 



0.29 0.85 



1.75 
2.08 
7.13 
0.03 
0.07 
0.15 
0.16 
0.62 



1.77±0.15 
2.04±0.25 



0.23±0.25 



Se 



BCS 
5.57 
5.53 
7.90 
0.01 
0.02 
0.07 
0.12 
2.78 



Q 

5.05 
5.00 
7.54 
0.04 
0.14 



S 
4.12 
3.63 
7.37 
0.08 
0.16 



exp 
4.4±0.2 
3.8±0.4 



0.27 0.39 
0.56 0.58 
3.55 5.66 



5.8±0.3 



2.77 
2.95 
7.76 
0.00 
0.01 
0.02 
0.03 
0.46 



2.76 2.28 
2.97 3.03 
7.12 7.06 
0.03 0.04 
0.09 0.09 
0.17 0.18 
0.31 0.27 
1.15 1.04 



2.08±0.15 
3.16±0.25 



0.84±0.25 



To further test the adequacy of the adjusted effec- 
tive mean field, the running sum of the contributions to 
the M^"^, the nuclear matrix element for the 2z//3/3-decay 
mode, is shown in Fig. [3] and compared with the available 
data. The quantity displayed is 



{f\\c 



■+ll4)(l^|l.t)(l,tJ|g 

Wm-(M,+M/)/2 



(12) 

where on the x-axis in Fig. [3]we use the excitation energy 
Eex in the intermediate nucleus ^^As instead of the 
the largest included eigenvalue of the QRPA equation of 
motion. 

We adopted the adjusted effective mean field energies 
shown in Figs. [l]and|2]for ''^Ge and ^^Se together with 
the SRQRPA method for the evaluation of the 2i/(3f3 nu- 
clear matrix clement. Again, the adjusted effective mean 
field describes much better the energies and amplitudes of 
the states for excitation energies below ~4 MeV where 
experimental data are available. We should point out, 
however, that the good agreement for the product of the 
two weak amplitudes (the numerator of Eq.( 12 1) does not 



mean that our calculation is free from the general prob- 
lem of QRPA calculations, namely that the (3~ strength 
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corresponding to the GT transition ''^Ge ^^As is too 
large while the P'^ strength ^^Se ^^As is too small, as 
stressed, e.g., in Ref. |T8]. Based on Fig. [sjone can con- 
clude that the reasonable agreement between the experi- 
mental value of Af^"^ and the value based on the few low- 
lying states (the so-called low-lying states dominance) 
appears to be accidental, at least in this case; if mea- 
surements could be extended to ~5-6 MeV and stopped 
there, that agreement would be lost, particularly when 
the adjusted energies are used. 




FIG. 3; The running sum of M^", Eq. The dashed hue 

corresponds to Refs. [3], the full line is the present result, 
and the dot-and-dashed line is the experimental result [17j . 
In the insert the first 4 MeV of excitation energy are shown 
in detail. 

While the chosen procedure is not ideal, since the ef- 
fective single-particle energies are chosen ad hoc, it is 
clearly an improvement when compared to our previous 
work [3|. First, the adopted method, SRQRPA, conserves 
the mean particle number in the correlated ground state, 
unlike the QRPA. Second, the occupancies of the valence 
orbits, in both protons and neutrons, agree now within 
errors with the experimental data [51 [6]. 



IV. RESULTS 

We are now in the position to ascertain to what degree 
the Oj^/3/3-decay nuclear matrix element M^" for the ^^Ge 
—> ^^Se transition change with the modification of the 
eS'ective mean field energies in combination with the ap- 
plication of the SRQRPA method. In Table |III| we show 
the sequence of the M'^'^ values, together with the fitted 
renormalization constants gpp as well as the calculated 
halflife for the nominal neutrino mass (mpp) — 50 meV. 
The Table entries were obtained using ^°Ca as a core, 
with the p,f and s,d,g shells (9 single particle orbits) 
for both neutrons and protons included, with unrenor- 
malized value = 1.25, CD-Bonn based G-matrix in- 
teraction, and with Jastrow-like function included 19J in 
order to include the effect of the short range correlations. 



(We will discuss the dependence on these assumptions 
below. Note that since we use gA = 1.25 everywhere in 
this work, we denote the OvPfi nuclear matrix element as 
M^" and do not use the notation M'^" = (5^/1.25)2M°'' 
of Refs. [3^.) 



TABLE III: The values of gpp, M^"" , and T^J^ for (m^/s) = 50 
meV in units of lO^^y, for different variants of the calculations. 
For further explanations, see text. 



variant 


9pp 




ri72(io-y) 


QRPA, WS mean field [3] 


0.849 


4.34 


7.0 


RQRPA, WS mean field [3] 


0.990 


3.81 


9.1 


QRPA , effective mean field 


0.903 


4.23 


7.4 


RQRPA, effective mean field 


1.170 


3.44 


11.1 


RQRPA, effective mean field 
and adjusted M^'^ denominators 


1.206 


3.18 


13.0 


SRQRPA, effective mean field 


1.125 


3.66 


9.8 


SRQRPA, effective mean field 
and adjusted M^" denominators 


1.184 


3.27 


12.3 



Th e las t, and thus our final entry for this variant in 
Table III M'^'^ = 3.27, was obtained when, in addition, 
the lowest energy denominator in the expression for A/^'', 



Eq.( 12 ), was shifted in such a way that it agreed with the 
known energy of the first l"*" state in ^^As. All other en- 
ergy denominators were then shifted by the same amount 
as the first one. (The entry in line 5 was obtained this 
way as well.) This procedure is commonly used when 
M'^'^ are evaluated, but was not employed previously in 
Refs. [3 . This value, M'^'^ — 3.27, can be compared now 
with the same quantity calculated by other methods. To 
make the comparison meaningful one has to keep in mind 
that Ai^"" contain R, the nuclear radius, as a factor. Un- 
fortunately, different authors use different conventions for 
R = To X A^^^. In our work we use rg = 1.1 fm. Pre- 
sumably the same is used in Ref. [10], even though the 
tq value is not explicitly quoted there, while in the shell 
model works [TT] Tq = 1.2 fm. 

With this correction included, the shell model value 
[11] of AfO" is 2.11, and the QRPA result of Ref. [ID] 
with the adjusted single-particle energies chosen there, is 
2.8. All these calculations use the same method, namely 
the Jastrow function [19, for the treatment of short range 
correlations. It, therefore, appears that the adjustment 
of mean field energies, in order to correctly reproduce the 
measured occupancies of the valence orbits [5] |5] , results 
in reduction of M'^'^ when QRPA or its generalizations 
are used. The discrepancy with respect to the shell model 
result jllj is reduced, in our work, to about half of its 
previous magnitude. 

As in Refs. [3] we would like to estimate the possi- 
ble range of the M'^'^ values taking into account changes 
corresponding to the variation in the number of single 
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particle states included and to different treatment of the 
short range repulsion. Accordingly, we repeated the M^" 
evaluation with 3 and 5 oscillator shells included, in ad- 
dition to the two oscillator shell result described above. 
The single particle energies in these additional shells were 
kept at their original Woods-Saxon potential values. For 
all these three variants we performed the calculations us- 
ing the Jastrow fuction [TF| for the treatment of short 
range correlations as well as the Unitary Correlation Op- 
erator Method (UCOM) j2D]. We kept the axial current 
coupling constant qa at its unrenormalized value 1.25 in 
all cases. The previous, Refs. [5], and the new results 
are compared in Table IV 



TABLE V: The calculated occupancies of valence neutron and 
proton orbits for ^'^Se and *^Kr. See the text for explanation 





«^Se 

prev. new 


«^Kr 

prev. new 


neutrons 

P 

/5/2 
39/2 


5.5 5.2 
5.4 5.2 
7.2 8.5 


5.4 4.8 
5.4 4.9 
6.1 7.3 


protons 
P 

/5/2 
59/2 


2.7 2.2 
3.0 3.2 
1.0 0.9 


3.2 2.8 
3.6 4.0 

1.3 1.2 



TABLE IV: The calculated M"" matrix elements for the ^"^Ge 
Qvl3j3 decay; the mean value and its range are shown for the 
two alternative treatments of the short range correlations. In 
column 2 are the previous values obtained with (R)QRPA 
method and with the Woods-Saxon potential single particle 
energies [3, and in column 3 are the values obtained with the 
SRQRPA method and the adjusted energies described above. 





prev. new 


Jastrow s.r.c. 
UCOM s.r.c. 


4.24(0.44) 3.49(0.23) 
5.19(0.54) 4.60(0.39) 



V. APPLICATION TO THE ^'^ SE Qv^p DECAY 

It is reasonable to expect that the modifications of the 
mean field energies, described above, that were relevant 
to the ''^Ge —> ""^Se decay, will also apply, at least ap- 
proximately, to the ®^Se *^Kr decay, since in that case 
there are just two more protons and four more neutrons 
compared to the A — case. Guided by such consid- 
erations we modified the Woods-Saxon potential ener- 
gies used previously in Refs. [3 for ^ = 82 for the pro- 
tons in both ^^Se and ^^Kr as for the proton energies for 
^^Se, and for the neutrons as in the ''^Ge that is closer to 
them in the number of neutrons. The resulting valence 
orbit occupancies are shown in Table [V] In columns 2 
and 4 are the rfif''^^ values evaluated with the standard 
Woods-Saxon potential and in columns 3 and 5 are the 
"n(p)^^^^ values evaluated with the adjusted energies. In 
the neutron system the present treatment predicts that 
the valence shells contain only ~19(17) neutrons com- 
pared with the naive expectation of 20(18) neutrons. 

Performing the same set of calculations as described 
previously for the ^^Ge decay, we conclude that the Oi^/3/3 
nuclear matrix element for the ^^Se decay is also re- 
duced from the previous (R)QRPA value of 3.76(0.40) 
for Jastrow s.r.c. and 4.59(0.39) for UCOM s.r.c. to 
the SRQRPA values (with modified mean field energies 
and shifted energy denominators for the 2vf3P decay) of 
3.50(0.24) for Jastrow s.r.c. and 4.54(0.40) for UCOM 
s.r.c, again reducing somewhat the difference between 



this value and the shell model result of 2.0 |TT] (adjusted 
for the different values of rg and with Jastrow s.r.c). 



VI. CONTRIBUTION OF INDIVIDUAL ORBITS 

TO M"" 

In the QRPA, RQRPA, and SRQRPA the M^" is writ- 
ten as the sum over the virtual intermediate states, la- 
beled by their angular momentum and parity J'^ and 
indices ki and kf: 

(-i)^"+^-'+^+-^ X 

J'^ .ki,kf ,J' pnp'n' 

(pil),p^);J\\ /(ri2)OK/>i2) II n(l),n'(2);^) x 

{o^\\[4£n']j\\rkf){rkf\rk,){rkMc;cM\ot) . 

(13) 

The operators Ok,K = Fermi (F), Gamow-Teller (GT), 
and Tensor (T), contain neutrino potentials and spin and 
isospin operators, and RPA energies Ey^ ' . The J7'^ la- 
bels angular momentum and parity of the pairs of neu- 
trons that are transformed into protons with the same 



The nucleon orbits are labeled in Eq.( 13 1 by p, p' , n, n' . 
We can isolate the contribution of, say, neutron or- 
bits n,n' by fixing these two labels, but performing 
the summation over all other indeces. The result- 
ing two-dimensional array f(n,n') obviously must obey 
E„^„'/(7T., n') = M^'^; the individual contributions can be 
positive or negative. It is interesting to visualize such 
contributions in order to see which orbits are important 
and which are not, and to gain a better understanding of 
the various physics effects affecting the M*^" values. 

We show a lego plot of such contributions to the M'^'^ 
for the ^^Ge decay, normalized to unity, in Figs. |4] and 
[5] The large positive contributions along the diago- 
nals, stemming dominantly, but not exclusively, from the 
J'^ ~ 0+ pairing part of M""^, contribute -1-2.97 when 
added together. The off-diagonal entries, related to the 
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'broken pairs' or higher seniority parts of Al'-"^, give - 
1.97 when added. The well known opposite tendencies of 
the pairing and broken pairs contributions is thus clearly 
visible. In addition, one can also see that the valence 
orbits .g9/2iP3/2! /5/21P1/2 contribute considerably more 
than the orbits further away from the Fermi level, even 
though the and 37/2 give nonnegligible contributions. 




FIG. 4: The contribution of individual neutron orbit pairs to 
M°". The entries are normalized so that their sum is unity. 




FIG. 5: The contribution of individual proton orbit pairs to 
M"". The entries are normalized so that their sum is unity. 



lump all the other orbits further removed from the Fermi 
level into one combination r (for remote). This allows 
us to reduce the dimension of the matrix and the cor- 
responding plot to 10x10, shown in Fig. [6] Again the 
entries are normalized so that their sum is unity, and the 
labels along the x and y axes are arranged in such a way 
that most of the negative entries are in the front (total, 
naturally, again -1.97) and most of the positive entries 
are near the far corner, in order to enhance visibility. 

In Fig. |6]the contribution of the r non- valence remote 
orbits is sizable, and for the negative entries, in fact, dom- 
inating. However, the positive and negative contributions 
from combinations that include the r orbits cancel each 
other to a large extent (positive contributions total 1.29 
and negative ones -1.45) so that the net effect on M^'^ of 
the remote orbits is only ^15%. 




FIG. 6: The contribution of the initial neutron orbits (combi- 
nation n,n') (along the y axis) plotted against the analogous 
combinations p,p' (along the x axis) of the final proton orbits. 
The entries are again normalized to unity. For notation along 
the axes, see text. 

However, in QRPA and its generalizations the inclusion 
of orbits of at least two oscillator shells, i.e. the set that 
obeys the full Ikeda sum rule, is essential. Without it, 
the description of the 2i//3/3 decay is impossible with a 
reasonable value of the effective particle-particle coupling 
constant gpp. 



In order to better visualize which combinations of neu- 
tron and proton orbits contribute one could, in principle, 
isolate in Eq.(13l the pieces corresponding to the com- 
bination n,n' of neutron orbits from which the neutrons 
disappear in the initial nucleus and plot them against 
the combination p,p' of proton orbits in which the pro- 
tons appear in the final nucleus. Such a plot, however, 
would be difficult to visualize since it would represent a 
81x81 matrix even with our minimal space of nine or- 
bits. Instead, we consider just the three valence orbits 
p (representing both pi/2 and P3/2), /5/2 and 59/2, and 



VII. CONCLUSION AND 
ACKNOWLEDGMENTS 

The occupancies of valence neutron and proton orbits, 
determined experimentally in Refs.[5j |6], represent im- 
portant constraints for nuclear models used in the evalu- 
ation of the 0i>(3P nuclear matrix element. In the present 
work we have modified the input mean field in such a 
way that the valence orbits in the model obey these con- 
straints. Within QRPA and its generalizations we found 
that it is important to also choose the variant of the ba- 
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sic method that makes such comparison meaningful by 
conserving the average particle number in the correlated 
ground state. When following this procedure, but oth- 
erwise keeping the same steps as in our previous evalua- 
tion of AfO'^ within QRPA, we find that for the "^^Ge 
^^Se transition the matrix element is smaller by ~25%, 
reducing the previously bothersome difference with the 
shell model prediction noticeably. Moreover, when we 
assume that analogous changes in the mean field should 
be applied also to the ®^Se —^ ®^Kr Oi//3/3 decay, that dif- 
fers from the ^^Ge decay by only two additional protons 
and four additional neutrons, we find similar reduction in 
M^"^ as well. Clearly, having the experimental orbit oc- 
cupancies available, and adjusting the input to fulfill the 



corresponding constraint, makes a difference. It would 
be very useful to have similar constraints available also 
in other systems, in particular for ^-^"Te and/or ^"^^Xe. 
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